


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1969 


A study of an adaptive aircraft control system 
in a self-organizing configuration. 


Richarde, Henry Morgan 


Monterey, California. U.S. Naval Postgraduate School 


http://ndl.handle.net/10945/13134 


Downloaded from NPS Archive: Calhoun 


: Calhoun is the Naval Postgraduate School's public access digital repository for 
/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 











DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 


United States 
Naval Postgraduate School 


A STUDY OF AN ADAPTIVE AIRCRAFT CONTROL 


SYSTEM IN A SELF-ORGANIZING CONFIGURATION 


by 


Henry Morgan Richarde, Jr. 





October 1969 


This docwnent has been approved for public re- 
Lease and sale; its distribution 16 unloncted. 


Library 
U.S. Naval Postgraduate School 


Monterey, California 


93940 








A Study of an Adaptive Aircraft Control System 
in a Self-Organizing Configuration 


by 


Henry Morgan Richarde, erage 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1962 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
October 1969 


ABSTRACT 


An adaptive control scheme may provide the best approach to 
the problem of accommodating automatic flight control systems to 
the variations of dynamic characteristics encountered over the flight 
envelope of the aircraft. The C*-Criterion for aircraft time re- 
sponse provides a basis for the design of sucha system. 

To provide a margin of safety, all control systems have 
redundant channels for emergency use. An adaptive control system 
can be designed to be self-organizing, and in such a configuration 
can provide its own failure monitoring, resulting ina simpler and 
more efficient system. 

This study shows that such a system is feasible, that the 
response of the system is within the limits set forth in the C*- 
Criterion, and that the self-organizing characteristics provide 


reliable operation over the whole range of flight operations. 
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I. INTROS Ge EION 


A self-organizing system may be defined as one which changes 
its basic structure as a function of its experiences and/or environ- 
ment, with the general aim being to evolve toward some desired 
state of output behavior or mode of operation. Sucha system can 
be described as having the ability to learn, and can therefore be 
designed to handle situations which rigidly programmed systems 
are not particularly adept in handling. A report by Gibson, Fu, 
etal., Ref. 1, gives a good introduction to, and provides back- 
ground information on, learning control systems. A further and 
more recent introduction to the subject may be obtained from a 
report by Mendel, Ref. Z. 

Almost all of our present aircraft use some type of air-data 
scheduling as the method of accommodating their automatic flight 
control systems to variations in the dynamic characteristics en- 
countered over the full flight envelope of the aircraft. Such 
systems rely on measurements of dynamic pressure or Mach 
number. Also, each system must be tailored for the particular 
aircraft in which it will be deployed, with resulting extensive design 
changes from system to system, and expensive flight testing for 
final adjustments to insure adequate performance. 

With the rapid increase in technology in the field of electronics, 


many new electronic instruments have appeared and proven to be 
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almost indispensable in the safe operation of aircraft. These in- 
struments, already available and in use, can provide the necessa ry 
measurements for an adaptive flight control system. Sucha system 
would depend only on measurements of dynamic performance to 
provide the necessary gain-programming. This should give im- 
proved system performance as compared to a programmed- gain 
system. Since the self-adaptive technique uses a model reference, 
it should be possible to make sucha system almost universally 
applicable, requiring only a change in the model anda few parameter 
adjustments for each different aircraft. 

The method of approach used in this study is to extend that 
used for the North American SIDAC controller developed by 
Shipley and his associates, Refs 3 and 4. This was modified by 
Rang, Ref. 5, from the observation that the equation for the 
handling qualities given by the C*-Criterion, which was developed 
by Tobie, Elliot, and Malcom, Ref. 6, is very similar to the basic 
short period equation of motion of the aircraft. A feedback system 
using variable gains and fixed values representative of the ranges 
of the aircraft coefficients over the flight envelope were used to 
meet the C*-Criterion requirements. The gain changing mechanism 
is found by a gradient technique. The system is then recast into 
one which should show self-organizing properties. The development 
of the following theory is by Shipley [Ref. 4]. It is of the basic form 


and shows the essential ideas. The extension to a self-organizing 


configuration is given in the Appendix. 
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Beginning with the short period perturbation equations for 
straight and level flight, which may be written as 
1) B= mye + Mee Mek + MS, 
(2) n = U, Ce - «) =e eas 
angle of attack is eliminated because it is difficult to measure 


accurately and its measurement is affected by gusts and air tur- 


bulence. Equations (1) and (2) may then be written as 


(3) 9g = Mag e ae re | - “,. se ae 5 | + Mz; | 6 - “| 
a 
ia Ww Ze E ‘ ei - me 


The coefficients are abbreviated as 


LB = She, = me af 
BR. eee ea eye - Ma fy, 5 
By = Mo - 2,Ma/zn « 


These coefficients are negative for all flight conditions. Their 


values are listed in Table I. Equation (3) then becomes 


| OI 


(5) -£6 +36 +n + 2 S§ ae 1, 


The C*-Criterion for aircraft response is given by 
(6) C*¥ = wn +4,60 + 46 
This criterion defines the gravity forces that a pilot is subjected 


to as the aircraft responds to control inputs. 
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If C* is required to be a multiple of the command input C,, then 
(7) eal = R Ce 
and the criterion can be written as 


(8) Lo + U8 + n ~ RO = O : 


Equations (5) and (8) are of the same form and are exact if their 


coefficients are made to be equal. Thus 


dal] DI 


Ela — 
B= y, kc = 


ay vy 


de 


To make the problem as simple as possible, the effects of 
actuators, of any feedback through the aircraft control system, and 
of the dynamics of the control system were neglected. This ideal 
system is diagrammed in Figure 1]. It was then assumed that the 
elevator command was a linear addition of the stick command and 


the feedback. Then 

(9) lla a 

where the feedback was a function of the gains, so that 
CO) -  SveC a s ‘ 


and the gains © ; A , and as were such that there was no 


systemerror. The system error was defined as 
(11) € = © -~- 460 - Bn - AC 


When (10) was substituted into (9), and the resultant equation and 


(5) were together substituted into (11), the error became 
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Ga, ¢ = |B, - 4, + BM, | @ + | 2. -B, + AS.| Nn 
+ [Ap - 4, + 4%, | 
In order to hold the error to zero, the coefficients of 9, YI , 


and C must be zero. Therefore 


™ ESPEN C2 Barks C2 Bs - As 
a a 
Values of these coefficients with 
(oe r = > ‘3, = = ? B, i - 20 


are listed in Table II. 


The feedback system was composed of three separate feed- 
back calculators, each set for a given flight condition. These 
feedbacks were multiplied by gains K, and combined to give the 


total system feedback. The method of calculating these gains was 


developed by Shipley, [Ref. 4]. 


Letting 
WIM to * i ia ee 
where 

b, = 4, 6, ee, 

n By, Se A, ec 9 
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and K, , K, , and K, are positive numbers. 


Then differentiation with respect to time gave 


1s) 


DM ee 7 le al il ae dil dit 
(4) $Y = Bla Bd + BE + t B, Ate | 
and choosing 

ANG = 

At . K, 6G 

Zins = Ka ys G 

IRF 

AX's = CG 

7 Kee GS 


resulted in the equation 


Ch 2 
(15) SF = ALS EG 
Since B. <O, G was chosen so that 
Ue oe 
a O, -.<¢ € <€, 


al Ck en ae 


? 
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which insured that ay was negative. Therefore V_ will 


always decrease toward zero with time and the gains tended toward 


an ideal value for each flight condition. Inthis system the gains 


are calculated from the equations 
Cals Ke f @ sgn © dt 
T= Kafin sgne at 
We Kee sgn € lf 
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I EXPERIMENTAL. PROCEDURES 


The system was set up as shown in Figures 2 and 3, using 
the COMCOR 5000 analog computer. The basic aircraft analog 
was tested by making plots of © and VN asa step input Se was 
applied to the system. The values obtained were compared with 
computed values of steady state amplitude, frequency, and damping 
ratio. Transfer functions were derived for 3 and with respect 
to a step input, and the final value theorem was applied to obtain 
the steady state ratios. 

After the basic aircraft analog was tested, the feedback system 
was incorporated. Each of the three feedback elements was chosen 
to represent a particular flight condition. The three flight con- 
ditions chosen were 0002, 0009, and 5020, where the first two 
numbers in each series represent the altitude in thousands of feet 
and the last two represent the Mach number intenths. By choosing 
these particular flight conditions the operating range of the air- 
craft was fairly well represented. Flight condition 0002 is the 
low altitude, low airspeed, landing condition. Flight condition 
0009 is a subsonic, low altitude condition where the aircraft has 
excellent response to control inputs and a high damping ratio. 
Flight condition 5020 is a high altitude, high speed condition where 
the control response is not very good and the damping ratio is 


low. 
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The three feedback elements were tested by Setting the value 
of K for the element being tested at 1.0, and the K's for the other 
two elements at zero. Since the element was set up to represent 
exactly the flight condition being considered, the values of K 
should not change and the error should be zero. This proved to 
be the case for each of the three basic flight conditions when K 
was set inas an initial value. However, when operating the system 
with no initial values set in, the K's reached a steady state value 
which was somewhat different from the combination of a 1.0 and 
two zeros. Each flight condition did reach a steady state in which 
the error signal remained in its dead band. For anexample, see 
Fee Sal 0 laid). 

After both the aircraft analog and the feedback analog were 
tested, the system was tested for self-organizing characteristics. 
For simulation purposes, the values of b , , and C were all 
taken from the aircraft analog. Ina real system, these values 
would come from separate sensors. Sensor failures were simu- 
lated by changing the value of N , A ‘io A in one of the 
feedback elements from its calculated value to a different value. 


In this simulation, the values for errors were zero and one. 
(11) Maas pq at Ni, 10 ewe, C 
Setting the value to zero simulated the loss of signal input from 


a sensor, while setting the value to one simulated an erroneous 


signal input. One other type of error which was simulated consisted 
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of using a steady input in place of one of the variable signals for 

© , nN, and C from the aircraft analog. This would simulate 

a sensor which was stuck and producing a constant output, regard- 
less of the actual conditions. This last type error was not tested 
extensively, since the proper values for such signals were not 


known. 
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IV. Resets 


The full system, for all flight conditions tested, drove its 
gains to values such that the system error remained in its zero 
range. The time required for this system adjustment varied with 
flight condition, frequency and amplitude of the input commands, 
and the number of commands given. The flight characteristics 
did not vary significantly during this transient time. (See Figures 
135514, ander} 

The self-organizing characteristics of the system appear 
to be adequate. Simulated failures, as outlined in the method of 
testing, were compensated by variations inthe gains. The plots 
in Figures 16 thru 21 show little change in aircraft flight character- 
istics. In some cases, the value of the error signal could not be 
made to remain in its dead band, but always quickly drove toward 
zero. For the three design flight conditions, any simulated 
failures outside of the design feedback element showed no effect 
on flight characteristics. This is as expected, since the design 
element should handle the flight condition by itself. When the 
failure was simulated in the design feedback element, there was 
usually a big change in the variable gains, anda slight change in 
performance. The flight characteristics were still satisfactory. 

When multiple failures of the zero input type were simulated, 


the system was still able to function adequately. There was some 
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deterioration in the amplitude of response and the time response 

characteristics. The system error usually did not stay in its zero 
range when command signals were given, but always went rapidly 

to zero. The system was able to handle combinations of up to six 
simulated failures of inputs, as long as at least one each of the OS ; 
yn , and C inputs remained. 

When multiple failures were simulated using erroneous input 
signals, system reaction was not as good as with the zero input 
failure simulation. The system was still able to reach a semi- 
steady state condition in which the error would drive to zero, but 
response was generally poorer. The poorest results were obtained 
when the error simulation involved two of the same quantities, 
such as two ©@S. In such cases, the system could not tell which 
input was the correct one, and the resulting performance plots 
were quite bad. Similar results were obtained by mixing the 
modes of failure together. It is considered that part of the problem 
encountered was the result of testing techniques, in which the 
inputs were increased by as much as a factor of seven to provide 
the erroneous signals. 

One other method was used to simulate failures. A constant 
signal was used to replace one of the inputs. The system was able 
to handle this type of error very well. This type of failure was 
not tested extensively since representative error values were not 


known, 


Za 


V. CONCLUSIONS 


A self-organizing adaptive control scheme can handle gain 
scheduling within a control system in sucha Fea that the system 
response is relatively invariant as required by the C*-Criterion. 
The ability to alter system gains as necessary to meet the response 
criterion gives the system better response than can be realized 
with a fixed gain schedule system. The system is also able to 
accommodate sensor failures, thus providing a built-in backup. 
Since the system uses a model reference, it should be easily 
adaptable to all types of aircraft by changing the model and the 


fixed feedback parameters to those of the particular aircraft. 
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The equations of motion of the system are 


6= £9+2, n+ B & 


(1) | wre 
Oe \ OE Se ae. 
6é= @ - £9 - £,N= AE 


Changing the configuration from one with a single feedback f , 
which has three varying gains iS . : ee , to one which 
has three feedbacks of the form of 4 : 

». - t l 
Ve 9 +Vn + SS = 


Y, O tien + Yo S 


F, 


\ 


\ 


(2) os 
. - Se_+ Gin evec 


fmt with fixed gains,.. In matrix form 


OE 
(3) Y= ies a 
ee 


q - é 


is constant. The three feedbacks are combined with varying 


gains. 
1 
+ K = K 
(4) > 
K 
3 | 


to give the final form of the input into the elevator channel 


to 


(5) ££ = KF, + KF, + KF, 


To facilitate the derivation, change to vector-matrix notation. 


Let 
9 a be /, < 
i= 2 ee joey | B= | By on oe le 
C Be Bs \, 


The set of equations from (I!) then yields 


Se 2 Mae et wows 4.0 x 
(6) } ; i 
e > (2 - 4 ae 


where the prime denotes the transpose of the matrix. 


Consider the function 
(7) Deva Ge 2s ee, YOu es Ae) 


where Q is a positive definite matrix to be chosen later. 
Differentiation gives 


(8) AV = (2’/-B'+ & v')Q A st 


Since V is to be driven to zero, whichin turn means € is 


driven to zero, the requirement is 
Jd V 
Toe 


Thus, choose 


ie 
ox 
iI 
OD 
i 
G) 
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= 
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so that 


Ae 
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Since the parameter Be is always negative, taking 


eo = Sgn € 
makes 
© = 
At a p; |< | 
so that 
AV 
Tt g£ O 
If 
K O O 
(10) ar’ = oY Tee oO 
O OP 


then the system is the same as given in the body of the thesis. 


Equation (5) combined with (2) is 
(11) i Om 

But also from equation (1) 

(12) fF = ¥@ x 


Therefore 


y¥ = Vk 


Equation (9) is then 


AY | 7) ~1 
f° Ne a ee 
or 
AK _ s\7l 
(13) a | COU) se 
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Now choose 


Ao 
(ave =). "Ca 
so that 
£, 
a4) of 2 @' le le¢ 
F, 
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- K. O O 
al ae— O kK, O 
O Q Ks 
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i! 
Kw 
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(15) aT UK ALG 
dKs ee eG 
dt 


These produce the variable gains used in the self-organizing 


Syicteni.. once 


O) =e Ge 6 ew 


ie \. wie non-singular and Q, is positive definite, then Q is 
positive definite, which is necessary for concluding that € goes 


to zeroas V _ goes to zero. 
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FLIGHT 

CONDITION 0002 0009 1504 3006 4509 5020 

ALTITUDE Sac 15,000 | 30,000 | 45,000 | 50,000 
cr a 





[oe 132 313 1.5 | Bee2 | 85.8 102 


8, | | 40 |e |e | ose | se 
a 
; 1.08 | 


TABLE I ATRCRAFT PARAMETERS 


(F4 AIRCRAFT Units in radians, feet, seconds) 
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TABLE ITI FEEDBACK CONSTANTS 
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FIGURE 4 


Plots for flight condition 0002, showing the response of 
the basic aircraft alone and with feedback. 
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Plots for flight condition 0009, showing the response of 


the basic aircraft alone and with feedback. 
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FIGURE 6 


Plots for flight condition 5020, showing the response 
of the basic aircraft alone and with feedback. 
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FIGURE 7 


Plots for flight condition 1504, showing the response 
of the basic aircraft alone and with feedback. 
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Plots for flight condition 3006, showing the response 
of the basic aircraft alone and with feedback. 
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FIGURE 9 


Plots for flight condition 4509, showing the response 
of the basic aircraft alone and with feedback. 
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FIGURE 10 


Plot for flight condition 0002, showing K's driving 
toward their steady state values. 
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Plot for flight condition 0009, showing K's driving 
toward their steady state values. 
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FIGURE 12 


Plot for flight condition 5020, showing K's driving 
toward their steady state values. 
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Plot for flight condition 0002, showing the change in response 
as the K's go from zero toward their steady state value. 
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FIGURE 14 


Plot for flight condition 0009, showing the change in response 
as the K's go from zero toward their steady state value. 
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Plot for flight condition 5020, showing the change in response 
as the K's go from zero toward their steady state value. 
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FIGURE 16 


a 


Plots for flight condition 0002, with the value of ” 
and ¥ , respectively, set to zero. 
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FIGURE 17 
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Plots for flight condition 0002, with the value of 


and a , respectively, set 
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FIGURE 18 


Plots for flight condition 0002, with the value of ~ ee ; 
and ” , respectively, set to zero. 
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FIGURE 19 


Plots for flight condition 0002, with the value of Y, , Vo , 
and Ya , respectively, set to one. 
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PIGURE. 20 


Plots for flight condition 0002, with the value of vr. , es ; 
and 0. , respectively, set to one. 
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ght condition 0002, with the value of Ns 
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